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We study the conductivity of graphene with a smooth but particle-hole-asymmetric disorder 
potential. Using perturbation theory for the weak-disorder regime and numerical calculations we 
investigate how the particle-hole asymmetry shifts the position of the minimal conductivity away 
from the Dirac point e = 0. We find that the conductivity minimum is shifted in opposite directions 
for weak and strong disorder. For large disorder strengths the conductivity minimum appears close 
to the doping level for which electron and hole doped regions (“puddles”) are equal in size. 
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I. INTRODUCTION 

The fascinating electronic transport properties of 
graphene have been the focus of intense experimental 
and theoretical investigation over the recent years.^^ll^ A 
particularly remarkable feature is the appearance of a 
finite minimum in the electrical conductivity as func¬ 
tion of the doping level when the graphene sample is 
tuned close to the Dirac point where electrons behave 
as ultra-relativistic particles.® At the Dirac point, the 
density of states and, hence, the charge-carrier concen¬ 
tration vanish in the absence of disorder. In that case, 
the minimum conductivity has been shown to originate 
from evanescent modes in the bulk and its value is found 
to be CTo = 4e^/7r/i.®ti^ 

Disorder that is smooth on the scale of the lattice con¬ 
stant, a condition that is approximately met if the disor¬ 
der originates from charged impurities at a finite distance 
from the graphene sheet, increases the value of such mini¬ 
mal conductivity. This at first sight counterintuitive phe¬ 
nomenon can be understood from the observation that 
a smooth disorder potential always creates a finite car¬ 
rier concentration, independent of its sign. Therefore, 
disorder turns a uniform graphene sheet with zero car¬ 
rier density into a landscape of regions (“puddles”) that 
are electron or hole doped, separated by a network of 
carrier-free junctions pidHlSl xhe effect on the conductiv¬ 
ity results from a competition between this increase of 
carrier density and a simultaneous increase of scattering 
events among the charged particles. Analytical and nu¬ 
merical calculations have shown that increasing disorder 
indeed leads to an increase of the conduct ivity.C^HSO] 

The existing fully quantum coherent theories of trans¬ 
port near the Dirac point use simplified models for the 
disorder potential, in which the potential fluctuations are 
Gaussian. Such Gaussian disorder models are necessar¬ 
ily statistically particle-hole symmetric. The purpose of 
the present article is to investigate the effect of an inher¬ 
ent statistical asymmetry between electron and hole-like 
puddles. We note that such asymmetry is relevant for 
experiments where, e.g., the charged impurities in the 
substrate may have a preferred charge. Besides the ex¬ 
perimental relevance, the investigation of the influence of 


particle-hole-asymmetric disorder on the conductivity of 
graphene furthermore reveals information about the na¬ 
ture of charge transport in graphene which cannot be in¬ 
ferred from studies with particle-hole-symmetric (Gaus¬ 
sian) disorder. Finally, it has been realized recently, that 
disorder that does not break certain symmetries on the 
average may lead to so-called “statistical topological in¬ 
sulator” phases, which differ qualitatively from phases 
in which symmetries are not preserved in a statistical 
sense.l^ 

With particle-hole-asymmetric disorder, it becomes in¬ 
teresting to not only look at the magnitude of the minimal 
conductivity but also at its position as a function of the 
(Fermi) energy. Different models used to qualitatively 
describe electric transport in graphene predict different 
values for the doping level at which the minimal con¬ 
ductivity occurs. For example, if electro nic transport is 
determined by percolation physics, ^^^*231 Qj^g expects the 
conductance minimum to appear at a value of the en¬ 
ergy Esize for which the total sizes of particle- and hole¬ 
like puddles are roughly equal. On the other hand, in a 
mean-field Boltzmann approach, in which the disorder is 
characterized by an effective carrier density,the con¬ 
ductivity minimum would be found at an energy Echarge 
where the mean carrier density is zero. Finally, if the elec¬ 
trons are delocalized over distances that are much larger 
than the correlation length of the disorder, the charge 
carriers would feel only an effective potential averaged 
over a large area in space and the conductivity minimum 
would be expected at the energy e = 0 for which charge 
carriers feel a potential landscape with zero mean. The 
same location of the conductivity minimum appears, if 
one requires that the absolute carrier density is a min¬ 
imum. While all three values Egize, Echarge and e = 0 
are equal for particle-hole-symmetric potentials, they are 
in general different if the disorder potential lacks such 
symmetry. 

The considerations that we present in this article con¬ 
sist of two parts. First, in Sec.|TT]we will use an analytical 
perturbative approach suited to describe the conductiv¬ 
ity for finite system sizes and weak disorder strengths. 
Here we extend an e arlier study by Schuessler, Ostrovsky, 
Gornyi, and Mirlin,E2l2ll that addressed ballistic trans- 
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FIG. 1: Schematic picture of the setup considered for the con¬ 
ductance calculation: A short and wide {W S> L) graphene 
strip is connected to two highly doped leads. The leads are 
modeled by graphene strips with potential u —>■ oo. A volt¬ 
age V is applied to the contacts and we are interested in the 
resulting current I = GV. 


port in graphene with Gaussian disorder. Although the 
perturbation theory is limited to finite system sizes only, 
we find a size-independent shift of the position conduc¬ 
tivity minimum for weak particle-hole-asymmetric disor¬ 
der. In the second part, we analyze disorder of arbitrary 
strength using a numerical method. Our main result is 
that upon introducing particle-hole asymmetry the loca¬ 
tion of the conductivity minimum shifts in opposite di¬ 
rections for weak and strong disorder. For weak disorder, 
the numerical results are consistent with the perturbative 
approach. For strong disorder, we find that the conduc¬ 
tivity minimum roughly appears when the Fermi energy 
is tuned such that electron and hole puddles have the 
same size. Our numerical findings are discussed in Sec. 
mil We conclude in Sec. HV) Some technical details are 
provided in the appendix. 


II. PERTURBATIVE APPROACH 


Our analysis of the weak-disorder regime builds on pre¬ 
vious work by Schuessler et on the conductivity of 
graphene in the presence of particle-hole-symmetric dis¬ 
order. In this section, we first introduce the matrix Green 
function formalism used in Ref. 1251 and then apply the 
formalism to the case of weak particle-hole-asymmetric 
disorder. 


A. Matrix Green Functions and Conductance 

We consider a rectangular sheet of graphene with di¬ 
mensions L X VF in the short-and-wide limit W ^ L, 
see Fig. The sample is attached to source and drain 
contacts at a: = 0 and x = L and exposed to a smooth 
disorder potential that does not couple the valleys. The 


Hamiltonian reads 

H = rip • (T -I- U(r) — u{x) (1) 

where v is the velocity of the electrons in graphene, cr = 
(o'x,o'y) is a vector of Pauli-matrices corresponding to 
the pseudospin degrees of freedom for the two triangular 
sublattices of graphene, V (r) is the disorder potential and 
u(x) is the potential offset that accounts for the highly 
doped leadJiil 


u(x) 


0, 0 < X < L 

oo, else. 


( 2 ) 


To calculate the conductivity of the graphene sample 
we use the “Matrix Green function formalism”. This 
method, originally developed by Nazaro\EII and adapted 
for graphene by Titov et alW^, has been successf ully ap- 
plied in a variety of transport studies for graphene 
Central object is the matrix Green function (MGF) 
G(r,r'), which is defined as 


G-^ 


e — H -\-iQ —VxS{x) sin | \ 
—VxS{x — L)sin^ e — H — iO y’ 


( 3 ) 


where e is the electronic energy, Vx is the operator for 
the ^-component of the velocity, and </> is an additional 
parameter that will be set to zero at the end of the 
calculation. The parameterization with sin ^ is conve¬ 
nient for the MGF of clean and undoped graphene. For 
graphene, the Hamiltonian H and the velocity operator 
Vx = V Ox have an additional pseudospin structure such 
that G(r, r') is a 4 x 4 matrix. We further note that 
e = 0 corresponds to the case where the Fermi level is 
tuned to the Dirac point. A finite energy describes a 
doped graphene sample. 

In the Matrix Green function formalism, the conduc¬ 
tance G follows from the MGF G through the relation 


8e2 

h dcj)^ ^=0 


H = Tr In G 


( 4 ) 


where the trace Tr extends over spatial coordinates as 
well as pseudospin and retarded-advanced space. We 
have included a factor of 4 for spin and valley degeneracy. 
The conductivity a is obtained from the conductance G 
using the relation 


a = 


hm 

W>L W 


( 5 ) 


The limit W ^ L ensures independence of the transverse 
boundary conditions. 

Reference gives a closed but lengthy expression for 
the matrix Green function Go in the absence of disor¬ 
der, from which one obtains the well-known result for 
the clean conductanc^44l 


_ 4e2 W 
Trh L 


( 6 ) 
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FIG. 2: Diagrams that contribute to the conductance cor¬ 
rection. Crosses represent energy contributions and dashed 
lines correspond to disorder averages. Diagram (a) yields the 
second order energy correction to the graphene sample’s con¬ 
ductance. Diagram (b) is needed to determine the second 
order correction in the disorder potential. The two points for 
which the disorder average is performed are only separated by 
a small distance of order of the disorder correlation length 
Finally, (c) generates the conductance correction of a disorder 
potential’s particle-hole asymmetry. Here, three points which 
are close to each other have to be considered for the third 
disorder averaged moment. The energy acts at r' which does 
not have to be close to the other coordinates. The resulting 
conductance contribution is presented in equation (211. 


Taking the clean and undoped case as a reference point, 
the effect of doping and/or disorder can be systematically 
included using perturbation theory. Hereto, we introduce 
V = I^(r) — e as a perturbation operator that accounts 
for both doping and disorder effects. Correspondingly we 
write 


n = no + 5n, fIo = TrlnGo, (7) 


where the trace tr is performed over the 4x4 matrix 
structure of the MGF. This correction is represented by 
the diagram displayed in Fig. Hh- The evaluation of 
the trace and the spatial integrals in Eq. (10) yield the 
conductance correctioiP^ 


SGe = Cl 



( 11 ) 


with the numerical prefactor ci ~ 0.101. The same re¬ 
sult was found from an alternative calculation based on 
scattering theory.EI 

We now turn to the effect of a Gaussian disorder po¬ 
tential V (r) that has a zero mean and is characterized by 
the correlator 


(y(r)E(r')) 


a{hvY - 

271^2 ® 



( 12 ) 


where a and ^ label the strength and the correlation 
length of the disorder, respectively. We use angular 
brackets (...) to indicate the disorder average. For the 
description of graphene in terms of a single Dirac cone, 
it is essential to keep the correlation length ^ finite. To 
leading order in a, one finds for the correction to D 


^ J drdr' (tr [E(r)Go(r, r')E(r')Go(r', r)] ) 

a(hi;)2 f -iel/ 

X tr [Go(r,r-f p)Go(r + p,r)] (13) 


so that Do gives the clean-limit conductance Go, while 
i5D contains the corrections that exist due to disorder 
and doping. Explicitly, one has 


Sn = Tr 


oo . 


( 8 ) 


The latter equation generates diagrams that consist of 
closed loops of clean MGFs. As shown in Refs. 1251771 for 
such structures, one can work with a simplified version 

of Go, 


see Fig. for the corresponding diagram. We focus 
on the limit ^ L for which the correlation length is 
much shorter than the sample size. In this limit, one can 


expand the square bracket in the last line of Eq. (13) for 
small p. Concentrating on the term proportional to (f)^, 
which is relevant for the conductance, one finds that the 
p dependence of the product of the two MGFs in Eq. (13) 


can be neglected. The integrations with respect to p and 
r can thus be performed separately and yielcP^ 


5G„ = —Go. 

ZTT 


(14) 


Go{x,y]x',0) 


1 f i cosh^ sinh ^ \ 

dhvL \ sinh ^ —i cosh ^ / 


1 


sin 


1 


'T^{x—x'-\-iy<7z 


Sin 


2L 


( 9 ) 


Before discussing the case of particle-hole-asymmetric 
disorder, we briefly review some of the results for the case 
of Gaussian disorder, that were obtained by Schuessler et 
^^p5] First, we look at the effect of a finite energy e inside 
the sample in absence of any disorder. Here, the first 
non-vanishing correction arises at second order in £, 


5Vts = — 
2 


J dr y dr'tr [Go(r,r')G'o(r',i')] , (10) 


This demonstrates that smooth disorder leads to an in¬ 
creased conductance at the Dirac point. (Reference [53 
further calculates the correction oc a^. This result will 
not be repeated here). 


B. Particle-Hole Asymmetric Potentials 

Having reviewed the concept of a matrix Green func¬ 
tion and the results for Gaussian disorder, we now ap¬ 
ply this method to disorder potentials that have an in¬ 
herent statistical particle-hole asymmetry. Since a non¬ 
vanishing mean of the disorder potential would only cor¬ 
respond to a redefinition of the energy £, the lowest 
non-trivial effect of a particle-hole-asymmetric potential 
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arises from its third disorder-averaged moment. To be 
specific, we assume a correlator of the type 


r'|2 + |r'-r"|^ + |r"-r 


{V(t)V{v')V{v")) = ' 

(15) 

where /3 and A quantify the strength and correlation 
length of the third-order disorder correlator, respectively. 
Again, we take a finite correlation length into account 
to avoid any problems that arise at short distances for 
Dirac-like particles. For the second moment of the disor¬ 
der potential we continue to use Eq. (12). 


The effect of such particle-hole-asymmetric potentials 
can be investigated using perturbation theory follow¬ 
ing essentially the same steps as for the particle-hole- 
symmetric potential. First, we note that precisely at the 
Dirac point (e = 0) there will be no effect linear in /3 
since any odd order perturbation in V gives a vanishing 
contribution due to the particle-hole symmetry of the un¬ 
perturbed problem. In the vicinity of the Dirac point, the 
correlator (Eq. will give rise to a term proportional 

to e/3 that leads to a shift of the minimal conductivity 
away from the Dirac point. Calculating this contribution 
to D requires the evaluation of the diagram shown in Fig. 

Srip =—e J drdr'dpdp'K{r,r', p, p') 

X {V{y p)V{v)V{y + p')). (16) 

The function K contains the trace over four Green func¬ 
tions 

K{r,r',p,p') = tr [Go(r; r-H p)Go(r-K p; r') 

xGo(i’', r + p')Go{r + p'; r)] . (17) 


two spaces can be carried out separately. We find that 

1 {y - v’f 


§^K(vy.p.p') 


1 




cos — x') — cosh f (y — y') 
1 


cos f (a: + x') — cosh f (y — y') 


(19) 


up to sub-leading corrections of order p In Eq. (19) 
we abbreviated 

/(».0i.y)= + ,20) 

sin 9 sin 9i (cos^ 9 -|- cos^ 9i sin^ 9) 


Although the r.h.s. of Eq. ( ^ diverges proportional to 
for small p, we note tl 


four-dimensional sphere contains a factor p^ such that 


;hat the integration over the 

the radial integral is well-defined. 

We proceed with the calculation of the conductance 
by combining Eqs. (|^, (15), (16) and (19). Translational 
invariance in the transverse direction implies that one in¬ 
tegration over a y-coordinate simply results in the width 
W of our sample. The final result reads 


6Gp — —C2f3^r^Go, 


( 21 ) 


where the factor C 2 « 0.00221 can be expressed in the 
form of dimensionless integrals, see Appendix for more 
details. 


Equation (211 is the main result of this Section. As the 


change (21) is linear in energy, together with Eq. (11) it 


implies a shift of the position of the minimal conductiv¬ 
ity away from e = 0. Combining the three perturbative 
corrections to the conductance of graphene in the weak- 
disorder limit close to the Dirac point, we find 


We note that a factor ^ in the expansion of Eq. (§ is 
canceled by the four possibilities where e can be placed 
for the perturbation V. 

In the limit X L the main contribution from the 
disorder correlator comes from values close to p = p' = 
0. To perform the integrations, we parameterize p and 
p' with four-dimensional spherical coordinates, 

Px = pcos0, 

Py = p sin 6* COS 01, 

p'^ = p sin 9 sin 9i cos p, 

p() = psin0sin0i sinp. (18) 


~ Trh L 


a (sL 


( 22 ) 


where the numerical constants Ci and C 2 are given below 
Eqs. 0 and ( |21[ ) , respectively. From this expression, we 
find that the minimal conductance occurs at the energy 


^min Cmin 


phv 


(23) 


where Cmin = C 2 / 2 ci « 0.011. 


III. NUMERICAL APPROACH 


To obtain the leading correction of the conductance in the 
limit A ^ L, it is then sufficient to extract the behavior of 
the function K for small values of p, where we can focus 
on the terms proportional to Here, it is useful that 
Eq. (§ for Go is written as a tensor product in retarded- 
advanced and pseudospin space, so that the traces for the 


We now present a numerical calculation of the con¬ 
ductivity in graphene in the presence of particle-hole- 
asymmetric random potentials. The numerical simula¬ 
tions can access the regime of weak disorder as well as 
the regime of strong disorder. 
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FIG. 3: (Color online) Parabolic fit of the numerically ob¬ 
tained conductivity data. The parameters are a — 0.16/27r, 
L/^V2 = 40 and W/L = 5. 


Our calculations are based on a method introduced by 
Bardarson et al^ In this approach, the sample is sliced 
into many segments of dimension 5x xW where the ex¬ 
tension 5x is chosen small enough such that its scattering 
matrix is determined perturbatively using the Born ap¬ 
proximation. After concatenation of the individual scat¬ 
tering matrices, one obtains the transmission matrix t for 
the entire sample. The conductance G is then obtained 
from the Landauer formula 

4e^ , 

G=—^ tilth]. (24) 

For more details on the numerical method we refer to 

Ref. [T71 

To numerically generate a particle-hole-asymmetric 
disorder potential we start from an auxiliary random 
function U{r) with Gaussian fluctuations, (G(r)) = 0 
and 


(C/(r)C/(r')) 


j~a {hvY - 

V ^ 2^2 ® 



(25) 


We then construct the disorder potential V (r) as 


V(r) = 



I a {hvh 

V ^ 2^2 


hv 


(26) 


One easily verifies that (R(r)) = 0, that its second 
moment obeys Eg . (12), and that its third moment is 
given by Eq. (15), with parameters /? = (8/3)V27ra3 
and A = Note that the potential V[y) is bounded 
from below, but not from above, so that it is manifestly 
particle-hole asymmetric. 

Using the disorder model (26) we numerically calcu¬ 
lated the conductance as a function of energy, for var¬ 
ious values of the disorder strength a. We perform an 
average over at least 500 disorder realizations. The as¬ 
pect ratio lU/L = 5 is chosen large enough to ensure that 


FIG. 4: (Golor online) System-size dependence of the energy 
value Emin ut which the conductivity minimum is found. Up¬ 
per panel: disorder strength a = 0.04/27r. Lower Panel: 
a = 16/27r. The aspect ratio W/L = 5 in both panels. 
The disorder average was performed for 500 disorder real¬ 
izations at L/^y/2 = 40 and for 200 (25) realizations at 
L/^y/2 = 60 (100). 


the conductance is insensitive on the particular boundary 
conditions in the transverse direction. We then extract 
the energy emin at which the minimal conductivity of our 
graphene sample occurs with a parabolic fit. An exam¬ 
ple of such a fit is shown in Fig. i for L/^-\/2 = 40 and 
a = 0.16/27r. We verified that our result for emin does 
not strongly depend on L, see Fig. for representative 
traces for small and large disorder strengths. 

Figure shows the location of the conductivity mini¬ 
mum for small and intermediate disorder strengths (a 
1). For weak disorder, the numerical results are consis¬ 
tent with the perturbation theory of the previous Section, 
see the inset of Fig. Remarkably, beyond the weak 
disorder regime, the dependence of Emin on the disorder 
strength is non-monotonous, and Emin becomes negative 
for a > 1. This is shown clearly in Fig.[^ where we show 
£min for the full range of disorder strenghts a considered 
in the numerical simulations. 

To put these findings in perspective, we consider the 
three “naive” scenarios for the location of Emin mentioned 
in the introduction: If percolation physics describes the 
conduction in graphene close to the Dirac point, the con¬ 
ductivity is minimal at the energy Esize for which the 
sizes of particle and hole-doped regions are equal. If an 
effective carrier density is the relevant variable that de¬ 
termines electric transport in graphene close to the Dirac 
point, such minimum would be found at an energy Echarge 
for which the total charge of the sample would be zero. 
Finally, the minimum could also occur for e = 0 for which 
the landscape e — U(r) is zero on average. This is the 
point at which the absolute carrier density has a mini¬ 
mum. It is also the relevant energy if electrons are delo¬ 
calized over distances larger than the correlation length 
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and would only be affected by an averaged potential. 

We now calculate the energies Ssize and echarge for the 
random potential defined by Eq. (26). For such calcu¬ 
lation we can replace the spatial averages by ensemble 
averages. As the potential U{r) is Gaussian correlated, 
its probability distribution reads 


q{U) 




(27) 


The quantity Esize is then determined by the condition 
that the probabilities that the potential V (r) < ggize and 
that V (r) > Esize are precisely equal or, equivalently, that 


IV(U)<es: 


dUg{U) = 


(28) 


Using Eq. (26) for V(U) and evaluating the integral gives 
the condition 


Erf 




(29) 


with the error function Erf (a;) = (2/y^) fjdte so 
that 


e,ize^-0.154^^. (30) 

In a similar way, one can calculate the energy £charge 
for which the total charge of the system is zero. For a 
given value of the disorder potential V, the local carrier 
concentration n(£) is given by 



a 


1/2 


FIG. 5: (Color online) Location Emin of the conductivity mini¬ 
mum for weak and intermediate particle-hole-asymmetric dis¬ 
order. The data points are the result of numerical calculations 
with the disorder potential chosen according to the model 
(261. The result (231 of the weak-disorder perturbation the¬ 
ory is indicated with a black line. The inset gives a log-log 
plot showing the agreement between the numerical results and 
the perturbation theory. 


remain resistive and are likely to dominate the resistance. 
For weak disorder none of these three naive limits de¬ 
scribes the numerical results well. This is not a large 
surprise, since transport in graphene in this regime is in¬ 
herently quantum mechanical and the semiclassical pic¬ 
tures underlying the three naive estimates are likely to 
fail here. 


n(£) = sgn(£ — V) 


(s - f^)^ 

7r(hvy 


(31) 


IV. CONCLUSION 


Using this relation, we can calculate £charge from the con¬ 
dition 



dt/u(£charge)/?(b/) — 0, 


(32) 


where the potential V is given by Eq. (26). Solving for 

^charge giveS 


£charge 0.156^^. (33) 

Note that £charge fiss the same dependence on the dis¬ 
order strength and correlation length as £size, but that 
its sign is opposite. Inspecting Fig. we conclude that 
for large disorder the position of the minimal conduc¬ 
tivity agrees reasonably well with the energy scale £size 
for which particle and hole puddles equal each other in 
size. We explain the applicability of the “percolation- 
motivated” picture at large disorder strength with the ob¬ 
servation that the charge density in the puddles becomes 
large for strong disorder, whereas the puddle boundaries 


In this work, we studied how particle-hole asymme¬ 
try of the disorder potential affects the conductivity of 
graphene. Since there are no a priori reasons why a realis¬ 
tic disorder potential should be statistically particle-hole 
symmetric, it is desirable that simplified models for the 
disorder do not have a higher symmetry. Moreover, the 
analysis of the energy at which the conductivity mini¬ 
mum occurs (relative to the Dirac point £ = 0) reveals 
information about the transport mechanisms of graphene 
close to the Dirac point, which cannot be found in stud¬ 
ies with Gaussian disorder, for which the conductivity 
minimum necessarily occurs at the Dirac point £ = 0. 

For a small disorder strength a, we calculated £min 
from perturbation theory. This result was verified by nu¬ 
merical simulations. Our most remarkable observation is 
that upon increasing the disorder strength the sign of the 
energy £niin at which the conductivity minimum occurs 
changes. At strong disorder, a 1, the magnitude and 
sign of £„iin are consistent with a percolation picture, in 
which the resistivity is dominated by the resistive junc¬ 
tions between otherwise highly conducting electron and 
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FIG. 6: (Color online) Location Smin of the conductivity 
minimum for the full range of disorder strengths. The data 
points are the results of our numerical calculations, the lines 
represent the naive estimates Ssize (solid, black) and Scharge 
(dashed, green). For large disorder strengths, the numerical 
data shows a good agreement between emin and the energy 
value corresponding to equal sizes of electron and hole pud¬ 
dles inside the sample (eaize). 


hole puddles. 

Realistic values for th e disorder strength a are in the 
range 1 ^ a ^ 4 IIII 33 ] unfortunately, is in the 

range of intermediate disorder strengths, for which per¬ 
turbation theory fails, but the strong-a asymptotics has 
not yet set in. 
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Appendix: Calculation of the numerical constant C2 


In this appendix, we show some details on the calculation of the integrals that lead to the numerical prefactor C 2 of 
Eq. (21 1 . As explained in the main text, we transform the integrals over the displacements p and p' to four-dimensional 


spherical coordinates, see Eq. (18). The Jacobian for this variable change is 

r27r 


/ OO ^00 pTT pTi pAT\ 

dpdp'[...] = / dp dO d 0 i / d(/ 7 p^ sin^ dsin 0 i[...]. 

-00 Jo Jo Jo Jo 


Furthermore, the disorder correlator is expressed in terms of the spherical coordinates as 

{y(y- -I- ^)l/(r)l/(r -1- 5')) = cosysin2esinei-4 sinysin^ esin26li]/A^ 


(34) 


(35) 


In order to obtain Eq. (21), one utilizes Eqs. (|^, (15), (16) and (19) from the main text. Translational invariance 
in the y-direction implies that one j/-integral results in the system width W. We rescale the remaining y-integration 
as well as the two ^-integrations with L and the p-integration with A to obtain Eq. (21) from the main text. The 
numerical constant is then given by 


C2 — TpTr 


(36) 


where 


= 


r dp rde Td^i r d^ P sm e ^ ^ [l- ^ cos ^ sm 2^ sm e, - 1 sin ^ sin^ g sin 2^,] ^ 

Jo Jo Jo Jo COS^ 9 + cos^ 6i sin^ 9 


Ir = 


IdTT^ 


dy / dxdx' 


y 


y 


cosh Try — cos Tr{x — x') cosh ny — cos 7r(a; -|- x') 


We first discuss the calculation of Ip. The integration over the variable p is straightforward and gives 


Ip= d9 


ddi 


dp sind 


cos (fi cos d -I- sin p sin 9 cos 9i 


(2 — cos p sin 29 sin 9i — sin p sin^ 9 sin 29i) (cos^ 9 + cos^ 9i sin^ 9) 


(37) 


(38) 



















Proceeding with the iy9-integration, one writes 


cos Lp cos 9 + sin ip sin 9 cos 9i 
cos^ 9 + cos^ 9i sin^ 9 


= Re 


p*<p 


cos 0 + i sin 0 cos 9i 


Using the identity 


c 2 tz 


, 1 - \/l - _ /32 

dp - -—- = 27r(Q; + ijd) - =. 

l-Q;C0S(p-/3sinV3 (a 2 +/ 32 )y^l-Q; 2 _y 32 ' 


for |q;|, |/3| < 1 /2, we then find that 


pTT p'K 

Xp=Tr d9 d^iTTsin^ 0sin0i 

Jo Jo 


where we abbreviated 


1 -Vl-g(^,^i) 
gi9,9i)^l-g{9,9,y 

g{9, 9i) = sin^ 9 sin^ 9i (cos^ 9 + sin^ 9 cos^ 9i) . 


(39) 


(40) 


(41) 


(42) 


This integrand is now a well-behaved function of 9 and 9i and the remaining integrals can be performed numerically, 
with the result 


Ip « 5.67862. 

For the calculation of 1^ we start with the ^-integration, where we use the identity 


dy 


y 


2 <I)(<I) - 7r)($ - 27r) 


cosh y — cos 4) 3 


sin4> 


(43) 


(44) 


which is valid for 0 < 4> < 27r, and can be continued to the entire real axis via the periodicity in 4>. This leaves us 
with the integral 

1 j ! ^\x — x'\{\x — x'\ — l){\x — x'\ — 2) {x + x'){x-\-x'— l){x + x'— 2) 

"" Jg ^ ^ I sinTrja; —x'l sin7r(x-|-cc') 

For the remaining integrations over x and x', it is customary to shift to the sum and the difference of these variables. 
Since each term depends only on one of these quantities, the integration over the other variable is straightforward. 
After some algebra, one obtains 



Ir 


2'k'^ Jq sinTTM 

35C(3) 93C(5) 

Stt® 27r® 

0.000388585, 


where is the Riemann-C function. 
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